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Abstract 

If dark matter (DM) is unstable, in order to be present today, its lifetime needs to 
be longer than the age of the Universe, %~4x I0 17 s. It is usually assumed that if 
DM decays it would do it with some strength through a radiative mode. In this case, 
very constraining limits can be obtained from observations of the diffuse gamma ray 
background. However, although reasonable, this is a model-dependent assumption. 
Here our only assumption is that DM decays into, at least, one Standard Model (SM) 
particle. Among these, neutrinos are the least detectable ones. Hence, if we assume 
that the only SM decay daughters are neutrinos, a limit on their flux from DM decays 
in the Milky Way sets a conservative, but stringent and model-independent bound 
on its lifetime. 

PACS: 95.35. +d, 14.60.St, 95.55.Vj 

1 Introduction 

It is more than seventy years since an unknown missing mass was first postulated in order 
to understand the motion of galaxies in clusters PQ. Since then, a lot more evidences 
at different scales have been found in favor of the existence of this non-luminous matter, 
from the measurements of the rotation curves of galaxies to the observation of the cosmic 
microwave background (CMB) (for reviews see eg. Refs. |2J). Thus the question now is not 
if the dark matter (DM) exists, but what is its nature. Many non-baryonic candidates have 
been proposed, from the lightest particle in supersymmetric models, to the lightest Kaluza- 
Klein excitation in models of extra-dimensions, to sterile neutrinos, light scalar particles, 
axions, particles from little Higgs models, etc. (see eg. Refs. [2 J for a comprehensive list). 

It might well happen that DM consists of different species, with different interactions 
with Standard Model (SM) particles and among themselves. It might also happen that 
DM was not in thermal equilibrium in the early Universe, i. e. it is not a thermal relic. 
Nevertheless, in order for it to be present today, there is a requirement: it needs to be 
stable, or at least to have a lifetime longer than the age of the Universe, tjj ~ 4 x 10 17 s. 
The possibility of unstable DM is not new and models with decaying DM have been con- 
sidered since long ago for different aims in astrophysics and cosmology [3J. Bounds on 
its radiative lifetime have been usually obtained from measurements of the diffuse gamma 
ray background [4]. However, although very stringent, they do not necessarily represent 
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constraints on the actual DM lifetime, but only set limits on its decay mode into pho- 
tons. On the other hand, by evaluating how the expansion rate of the Universe is affected, 
bounds on the DM lifetime can be set using CMB data [6j. In this more general case, it 
is shown that already the first year of WMAP observations [5] constrains DM lifetime to 
be longer than tcmb = 1-6 x 10 18 s at 2cr confidence level (CL), with the only assumption 
that the decay is into relativistic particles [Bj. This bound is rather robust but difficult to 
improve by further observations, for the DM decay affects CMB only at large scales, for 
which errors are limited by cosmic variance [UJ . Nevertheless, a recent study [7] , which also 
takes into account type la supernovae (SN) data, improves this limit by about an order of 
magnitude, tqmb+sn = 2.2 x 10 19 s at 2cr CL. 

In this letter we obtain a lower bound on the DM lifetime, which is much more restrictive 
(several orders of magnitude) than that set by CMB and SN observations and model- 
independent, unlike that obtained from observations of the diffuse gamma ray background. 

Among the stable SM particles, neutrinos are the least detectable ones. Therefore, if we 
assume that the only SM products from the DM decay are neutrinos, a limit on their flux, 
conservatively and in a model-independent way, sets a lower bound on the DM lifetime. 
This is the most conservative assumption from the detection point of view, that is, the 
worst possible case. Any other decay channel (into at least on SM particle) would produce 
photons and hence would give rise to a much more stringent limit. Let us stress that this 
is not an assumption about a particular and realistic case. On the other hand, for the 
reasons just stated, it is valid for any generic model with unstable DM, which decays at 
least into one SM particle. Hence, the bound so obtained is a bound on the lifetime of the 
DM particle and not only on its partial lifetime due to the decay channel into neutrinos. 
Thus, following a similar approach to that of Refs. [HI El CHj, we consider this case and 
evaluate the potential neutrino flux from DM decay in the whole Milky Way, which we 
compare with the relevant backgrounds for detection: mainly the well known and measured 
atmospheric neutrino flux, which spans over about seven decades in energy. 



2 Neutrino Fluxes from DM Decay in the Milky Way 

In what follows we only study DM decays in the Milky Way and do not consider the diffuse 
signal from cosmic decays. In general, the latter is likely to be smaller than, or at most of 
the same order of, the former. 

If DM has a lifetime longer than the age of the Universe, r x > tu, the differential 
neutrino (plus antineutrino) flux per flavor from DM decay in a cone of half-angle if) around 
the galactic center, covering a field of view AQ = 2tt (1 — cosip), is given by 

d® = An j ^ R BcPo 1 dN 
dE u Atx m x T x 3 dE v * 

where m x is the DM mass, i? sc = 8.5 kpc is the solar radius circle, po — 0.3 GeV cm -3 
is a normalizing DM density, which is equal to the commonly quoted DM density at i? sc , 
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and J7aq is the average in the field of view (around the galactic center) of the line of sight 
integration of the DM density, which is given by 



2 7T 1 Z" 1 ['max 

Jaq = -t^s / / p{r)dld(cost/)'), 

Zaij rC sc Pq Jcosip Jo 



(2) 



where r = J R 2 C — 2Z_R SC cosip' + I 2 , and the upper limit of integration i 



is 



Z max = \j {Rl alo - sin 2 ijjR 2 c ) + R sc cos i/> . 



(3) 



This integral barely depends on the size of the halo -Rhaio, as long as it is larger than few 
tens of kpc, for the contribution at large scales is negligible. 

The neutrino (plus antineutrino) spectrum per flavor is given by dN/dE„. If DM is the 
lightest particle of the new sector, which is introduced to render a more complete theory, 
then it can only decay into SM particles. Hence, the most conservative assumption is that 
it decays into neutrino-antineutrino pairs, x ~ > y v. In this case dN/dE u = 2 5(E u — m x /2). 
However, if the lightest particle of the new sector (% L ) is stable, but the next-to-lightest 
particle (xnl) is long-lived, the latter could also constitute part of the DM and decay into 
the former plus one or more SM particles, which we conservatively assume to be neutrinos. 
Commonly, in this class of models, these two new particles are almost degenerate in mass, 
and thus the total energy of the produced neutrinos is approximately equal to the difference 
of their masses (AM). For two-body decays, xnl - ¥ Xl + v, dN/dE u = 6(E U — AM), 
whereas for three-body decays, xnl —> Xl + v+v, the neutrino (and antineutrino) spectrum 
is continuous with a maximum energy equal to AM. In what follows we shall consider 
for concreteness (and for comparison with the CMB bounds) the case of DM decay into 
neutrino-antineutrino pairs and obtain bounds on r x as a function of m x . From this 
calculation, it is straightforward to obtain a bound on the combination function 
of the neutrino energy, which for the second two-body decay case is equal to AM. There 
are two main points to take into account. Whereas in the first case there is a produced 
neutrino-antineutrino pair, in the second there is only one final neutrino (or antineutrino). 
On the other hand, in the second case only half of the DM decays (the next-to-lightest 
particle of the new sector), for the lightest particle of the new sector is assumed to be 
stable, and it also contributes (usually at comparable level) to the DM. This implies an 
overall factor of 4. Finally, although a detailed analysis for three-body decays is model- 
dependent, in general this case would give bounds, for a given neutrino energy, of the same 
order of magnitude of those for the two-body decay case. 

In Eq. ([1]), the factor of 1/3 comes from the assumption that the decay branching ratio 
is the same for the three neutrino flavors. Let us note that this is not a very restrictive 
assumption, for even in the case DM decays predominantly into one flavor, there is a 
guaranteed flux of neutrinos in all flavors thanks to the averaged neutrino oscillations 
between the source and the detector. Hence, although different initial flavor ratios would 
give rise to different flavor ratios at detection, the small differences affect little our results 
and for simplicity herein we consider equal decay into all flavors. 
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On the other hand, with our definition of J ma, all the astrophysical uncertainties in 
the calculation of the neutrino flux from DM decays are encoded in J^aq- They come from 
our lack of knowledge of the exact DM density p(r). As a matter of fact, the formation of 
large scale structure is successfully predicted by detailed N-body simulations which show 
that cold DM clusters hierarchically in halos. The simulated DM profile in a galaxy like 
the Milky Way, assuming a spherically symmetric matter density with isotropic velocity 
dispersion, can be parametrized as 



where p sc is the DM density at R SC} r s is the scale radius, 7 is the inner cusp index, (3 is 
the slope as r — > 00 and a determines the exact shape of the profile around r s . 

The three commonly used DM density profiles we consider []^[T3J[TT] (see also Ref. [T4"] ) 
tend to agree at large scales, but uncertainties can be significant in their inner parts. 
However, for a large field of view, these uncertainties are much less relevant and do not 
affect significantly the calculation of the neutrino flux from DM decay In addition, and 
unlike the case of DM annihilations, the neutrino flux depends on the line of sight integral 
of the DM density and not of its square, which reduces considerably the effect of the inner 
cusp uncertainty. 

As we will see below, and following a similar approach as in Ref. [9], we are interested 
in signals corresponding to different components of the halo: the full-sky signal and the 
signal from a 30° half-angle cone around the galactic center. Whereas for the former, the 
value of the average of the line of sight integration of the DM density, J^iso, for the three 
considered profiles, can vary at the very most from 1.3 to 8.1, for the latter, the value of 
J730 might be anything from 3.9 to 24. These limiting cases are obtained from the range 
of values for p sc [T5] which satisfy present constraints from the allowed range for the local 
rotational velocity [16] , the amount of flatness of the rotational curve of the Milky Way 
and the maximal amount of its non-halo components [T7J . For the usually quoted value of 
p sc , for each of the profiles, (p sc , Ji 80 , J30) = (0.27 GeV/cm 3 , 1.9, 6.5) [11], (0.30 GeV/cm 3 , 
2.0, 6.1) [12] and (0.37 GeV/cm 3 , 2.2, 5.5) [13]. Thus, uncertainties in the halo profile have 
fairly small effects on our final results. Here we consider the simulation by Navarro, Frenk 
and White (NFW) [12] as our canonical profile. From the limiting cases just discussed, 
this implies that in the worst scenarios we could be overestimating (underestimating) the 
neutrino flux by a factor of about 1.5 (3.9). 

3 Neutrino Bounds on the DM lifetime 

For E v ^ 50-60 MeV, the main source of background for a possible neutrino signal from DM 
decays is the flux of atmospheric neutrinos, which has been measured in a number of detec- 
tors up to energies of ~ 100 TeV [H] • Its spectrum is also well understood theoretically and 
different calculations using different interactions models agree within 20-30% [T9" | |2"01 |2"T] . 
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Thus, in order to obtain a bound on the DM lifetime we need to compare these two fluxes, 
and in particular we consider the + spectra calculated with FLUKA |21j. 

Assuming two-body DM decays into neutrino-antineutrino pairs, we first obtain a gen- 
eral lower bound for m x ~ 100 MeV-200 TeV, by comparing the (z/ M + V^) neutrino flux 
from DM decays in the halo with the corresponding atmospheric neutrino flux in an energy 
bin of width Alog 10 E u = 0.3 around E v = m x /2. For each value of m x , the limit on r x is 
obtained by setting its value so that the neutrino flux from DM decays in the Milky Way 
equals the atmospheric neutrino spectrum integrated in the chosen energy bin. The reason 
for choosing this energy bin is two-fold: on one side, the neutrino signal is sharply peaked 
around a neutrino energy equal to half of the DM mass and this choice is within the exper- 
imental limits of neutrino detectors, and on the other side, for the sake of comparison, we 
follow the approach of Ref. 0. Nevertheless, following Ref. [TO] , a more detailed analysis 
is performed below for m x ~ 30-200 MeV. 

The most conservative bound is obtained by using the full-sky signal, and this is shown 
in Fig. [I] where the dark area represents the excluded region. However, a better limit can 
be obtained by using angular information. This is mainly limited by the kinematics of 
the interaction. In general, neutrino detectors are only able to detect the produced lepton 
and its relative direction with respect to the incoming neutrino depends on the neutrino 
energy as A9 ~ 30° x GeV/E u . As in Ref. [9] and being conservative, we consider a 

field of view with a half-angle cone of 30° (30° x ^10 GeV/E u ) for neutrinos with energies 
above (below) 10 GeV. This limit is shown in Fig. [T]by the dashed line (light area), which 
improves upon the previous case by about a factor of three for m x > 10 GeV. 

As anticipated, it is expected that a more detailed analysis, making a more careful use of 
the directional as well as energy information for a given detector, will improve these results. 
Note for instance that for energies ~ 1-100 GeV neutrino oscillations would give rise to 
a zenith-dependent background, roughly speaking a factor of two larger for downgoing 
neutrinos as compared to the upgoing flux, whereas we expect a nearly flat background 
for other energies for which oscillations do not take place. On the other hand, the signal 
from DM decays in the halo is expected to change by a factor of ~ 7 for a half-angle cone 
of 30° around the galactic center as compared to the signal within the same field of view 
but coming from the opposite direction. Thus, making use of the directional information 
would certainly render more stringent bounds. Nevertheless, and although a detailed and 
detector-dependent analysis is beyond the scope of this letter, we show how such a more 
careful treatment of the energy resolution and backgrounds can substantially improve these 
limits. For this and also extending the bounds to lower DM masses, we consider the low 
energy window below ~ 100 MeV (i. e. for m x ^ 200 MeV) and perform an analogous 
analysis to that in Ref. [TO] . 

In this energy range the best data comes from the search for the diffuse supernova 
background by the Super-Kamiokande (SK) detector which has looked at positrons (via the 
inverse beta-decay reaction, V e +p —> e + +n) in the energy interval 18 MeV-82 MeV [22]. As 
for these energies there is no direction information, we consider the V e signal coming from 
the whole sky. In this search, the two main sources of background are the atmospheric v e 
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m % [GeV] 

Figure 1: Bounds on the DM lifetime for a wide range of DM masses obtained using 
different approaches: full-sky signal (dark area), angular signal (light area) and 90% CL 
limit using SK data at low energies [22] (hatched area). Results are obtained for a NFW 
profile and assuming two-body decays into relativistic particles (see text). Also shown the 
bound obtained from CMB observations [6] and CMB plus SN data [7] (both at 2a CL) 
and the line r x = t\j. 



and z/ e flux and the Michel electrons and positrons from the decays of sub-threshold muons. 
Below 18 MeV, muon-induced spallation products are the dominant background, and below 
~ 10 MeV, the signal would be buried below the reactor antineutrino background. 

Although for E v ^ 80 MeV the dominant interaction is the inverse beta-decay reaction 
(with free protons), the interactions of neutrinos (and antineutrinos) with the oxygen nuclei 
contribute significantly and must be considered. For our analysis we have included both the 
interactions of V e with free protons and the interactions of v e and V e with bound nucleons, 
by considering, in the latter case, a relativistic Fermi gas model [23] with a Fermi surface 
momentum of 225 MeV and a binding energy of 27 MeV. We then compare the shape of 
the background spectrum to that of the signal and perform a x 2 analysis so that we can 
extract the limit on the DM lifetime in an analogous way as it was done to obtain an upper 
bound on the annihilation cross section for the case of DM annihilation in Ref . [TU] , where 
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we refer the reader for all the details of this analysis 0. The 90% CL limit is shown in 
Fig. [T] by the hatched area and it clearly improves (and extends to lower masses) by up to 
an order of magnitude upon the general and very conservative bound obtained with the 
simple analysis described above. 

Finally, let us note that in principle, if the DM mass is not known, DM annihilation and 
DM decay in the halo might have the same signatures. However, whereas the decay signal 
depends linearly on the DM halo density, the annihilation signal depends on its square. 
Hence, in case of a positive signal, directional information is crucial to distinguish between 
these two possibilities. 

4 Conclusions 

In this letter we have obtained a general bound on the DM lifetime, which is several orders of 
magnitude more stringent than previous limits [El [7]. In order to do so, we have considered 
that the only SM products from DM decays are neutrinos, which are the least detectable 
particles of the SM. Thus, regardless of how likely this is, by making this assumption we 
can obtain a conservative but model-independent bound on the DM lifetime. To do so 
we have considered the flux of neutrinos coming from DM decays in the Milky Way for 
an energy interval from ~ 50 MeV to ~ 100 TeV and have compared it to the dominant 
background, the well known and measured atmospheric neutrino flux. For concreteness 
we only show results for two-body DM decays into relativistic SM particles, although it 
is straightforward to generalize this result to other two-body decays. On the other hand, 
the model-dependent case of three-body decays should render limits of the same order of 
magnitude. We have obtained a general bound by considering the signal from the whole 
Milky Way and imposing that it has to be at most equal to the background in a given 
energy interval. We have also shown how this crude, but already very stringent limit, can 
be substantially improved by more detailed analysis which make careful use of the angular 
and energy resolution of the detectors, as well as of backgrounds. In this way, following 
the analysis of Ref. [10J, we have obtained the 90% CL lower bound on the DM lifetime 
for m x ~ 30-200 MeV, which is an order of magnitude more stringent. 

In summary, we have shown that neutrinos can be used to set very stringent and model- 
independent bounds on the DM lifetime, with the only assumption that if DM is unstable, 
it decays into at least one SM particle. As our main result, we have improved by several 
orders of magnitude upon previous limits [HI [7] . 
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